Measuring the Speed of Dark: Detecting Dark Energy Perturbations 
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The nature of dark energy can be probed not only through its equation of state, but also through 
its microphysics, characterized by the sound speed of perturbations to the dark energy density and 
pressure. As the sound speed drops below the speed of light, dark energy inhomogeneities increase, 
affecting both CMB and matter power spectra. We show that current data can put no significant 
constraints on the value of the sound speed when dark energy is purely a recent phenomenon, but 
can begin to show more interesting results for early dark energy models. For example, the best fit 
model for current data has a slight preference for dynamics (w(a) ^ — 1), degrees of freedom distinct 
from quintessence (c s 7^ 1), and early presence of dark energy (fide (a < 1) / 0). Future data may 
open a new window on dark energy by measuring its spatial as well as time variation. 



I. INTRODUCTION 

Although dark energy dominates the energy density 
of the universe and drives the accelerating cosmic ex- 
pansion, we know remarkably little about it. Over the 
course of the past decade, cosmologists have devoted con- 
siderable effort to devising new and sharpening known 
methods for determining the equation of state of dark 
energy. The equation of state, defined as the pressure to 
energy density ratio, is generally a time dependent func- 
tion and fully specifies the temporal evolution of dark 
energy density. The dark energy density in turn (along 
with the matter density) determines the expansion rate of 
the universe, as well as geometrical measures (distances 
and volumes). 

The equation of state w(z) does not, however, tell us 
about the microphysics of dark energy, nor does it de- 
scribe all of the cosmological signatures. For example, 
even a perfectly measured w(z) does not tell us whether 
dark energy arises from a canonical, minimally coupled 
scalar field, a more complicated fluid description, or mod- 
ification of gravitational theory on large scales. The prop- 
erties of the perturbations to the dark energy, which must 
exist unless it is simply a cosmological constant or only 
an effective field, do carry such extra information. 

Perturbations to the energy density and pressure can 



be described through the sound speed, 



5p/5p. 



The sound speed carries information about the inter- 
nal degrees of freedom: for example, rolling scalar fields 
(quintessence) necessarily have sound speed equal to the 
speed of light, c s = 1. Detection of a sound speed distinct 
from the speed of light would indicate further degrees be- 
yond a canonical, minimally coupled scalar field. 

A low sound speed enhances the spatial variations of 
the dark energy, giving inhomogeneities or clustering. 
Heuristically, the sound speed determines the sound hori- 
zon of the fluid, l s = c s /H, where H is the Hubble 
scale. On scales below this sound horizon, the fluid is 
smooth; on scales above l s , the fluid can cluster. Since for 
quintessence c s — 1, the sound horizon equals the cosmic 



horizon size and there are essentially no observable inho- 
mogeneities. However, if the sound speed is smaller, then 
dark energy perturbations may be detectable on corre- 
spondingly more observable (though typically still large) 
scales. These perturbations act in turn as a source for 
the gravitational potential, and affect the propagation 
of photons. For example, clustering dark energy influ- 
ences the growth of density fluctuations in the matter, 
and large scale structure, and an evolving gravitational 
potential generates the Integrated Sachs- Wolfe (ISW) ef- 
fect [l[ in the cosmic microwave background. The obser- 
vational signatures of these effects offer a way of probing 
the dark energy inhomogeneity and sound speed. 

In this paper we study the signatures of the sound 
speed of dark energy. We revisit and extend previous 
studies of dark energy clustering @-[2lj|, clarifying and 
quantifying the physical effects caused by the nonstan- 
dard values for the speed of sound. We then study models 
where the dark energy density was non-negligible at early 
times, which offer much better prospects for observable 
c s signatures than the fiducial near-ACDM case. Finally, 
using current cosmological data, we constrain the speed 
of sound jointly with 7-8 other standard cosmological pa- 
rameters. 

This paper is organized as follows. In Sec. |H] we de- 
scribe dark energy perturbations and the physical influ- 
ence of the sound speed and equation of state, deriving 
the dark energy density power spectrum. Section |TTT] de- 
scribes the dark energy models we consider, and Sec. IIVI 
treats the impact of dark energy inhomogeneity on the 
CMB, matter power spectrum, and their crosscorrelation. 
We consider models with both constant and time varying 
equation of state and sound speed in Sec. El and present 
constraints from current data. 



II. DARK ENERGY PERTURBATIONS 

We briefly review the growth of density perturbations, 
in both the matter and dark energy, focusing on the role 
of the sound speed. See 0, [23| for more details. To de- 
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rive the influence of the sound speed on dark energy inho- 
mogeneity, and dark energy perturbations on the matter 
distribution, we must solve the perturbed Einstein equa- 
tions for the density perturbations 5pi, pressure pertur- 
bations Spi, and velocity (divergence) perturbations 9i. 
We do not consider an anisotropic stress. 

In the conformal Newtonian gauge, the perturbed 
Friedmann-Robcrtson- Walker metric takes the form 



ds 2 = a(r) 2 [-(1 + 2tP)dr 2 + (1 - 2<f>)dr : 



(1) 



where a is the scale factor, r is the conformal time, r 
represents the three spatial coordinates, and ip and <j> are 
the metric potentials. Conservation of the stress-energy 
tensor (T^„ = 0) of a perfect fluid gives the following 
equations in Fourier space (see, e.g., [23() from the time- 
time and space-space parts: 



1 + w \op J 1 + w 



(2) 



= -ft(l-3iu)0- 



Sp 



+ ^k 2 — \-k 2 ip, 



1 + w Sp 1 + w 



where k is the wavevector, dots are derivatives with re- 
spect to conformal time, ft = a/a is the conformal Hub- 
ble parameter, 5 = 5p/p is the density perturbation, 
(p + p) 9 = ik°8T j is the velocity perturbation, and 
w = pi p is the equation of state. These equations hold 
for each individual component, i.e. matter or dark energy. 

We define the effective (or rest frame) sound speed c 8 
through (see, e.g., [24j |) 



S -P=clS + 3H(l + w)(c 2 s -c 2 a )^ 
p k 

where the adiabatic sound speed squared is 

2 _ P 1 w 

Ca = p~ = w ~mT+^i- 

In terms of c s , Eqs. and ^ read 



m{w - c 2 ) 



1 + w 



1 + w 



(3) 



(4) 



(5) 



-[k 2 + m 2 (c 2 s -c 2 a )} -2+30 



k 2 



(3c 2_l)^ + c 2_A_ + ^ (6) 



One can readily see that the source term in a 6 equa- 
tion will have a negative term involving c 2 k 2 from 9 (take 
the derivative of Eq. [5] and substitute in Eq. [6]), indicat- 
ing that growth is suppressed on small scales, k > ft/c s . 
However, perturbations will exist in the dark energy den- 
sity even for c s = 1, albeit at a very low level within the 
Hubble scale k > ft. As c s drops below unity, the sup- 
pression is itself suppressed and inhomogeneities in the 
dark energy can be sustained. All such perturbations 
will vanish though as 1 + w — > 0, regardless of c 2 . In 



the combination of Eqs. §5§ and ((6j) into a single second 
order equation for 5, the terms involving the metric in 
this equation are all proportional to 1 + w (or derivatives 
thereof) so that in the limit 1 + w —> the perturba- 
tions decouple from the metric and do not experience a 
gravitational force leading to growth. 

The dark energy perturbations affect the metric per- 
turbations, and thus the perturbations in the matter, 
through the Poisson equation 



k 2 4> = -4nGa 2 ^ Pi ( 5, + 3ft (1 



where the sum runs over all components. For a perfect 
fluid, there is no anisotropic stress so ijj = <j>. 

Therefore we expect the density power spectrum to 
be affected by the dark energy sound speed in distinct 
ways on different scales. On supcrhorizon scales, k < ft, 
the density power spectrum becomes independent of the 
dark energy sound speed. Here the perturbations are de- 
termined by the curvature fluctuation (25l . [26j . Between 
the Hubble scale and the sound horizon, ft < k < ft/c s , 
a sound speed c s < 1 will enhance the density inhomo- 
geneities (modulo gauge dependence around the Hubble 
scale). Finally, on smaller scales, k > ft/c s , inhomo- 
geneity growth is always suppressed and the sound speed 
becomes irrelevant. We illustrate these behaviors in Fig- 
ure Q] (All power spectra in this paper are for linear the- 
ory and shown at a = 1, and are calculated using CAMB 
[27j and CMBeasy [H, H].) Note that the strength of 
the deviation from the c s = 1 behavior is a steep function 
of c s for c s ps 0.1. 



III. DARK ENERGY MODELS 

We study three classes of dark energy models to elu- 
cidate the role of sound speed and 1 + w, from early to 
late times. 

1) Constant w models. We begin with the sim- 
plest model of dark energy with sound speed different 
from the speed of light: a constant equation of state w 
and a constant sound speed c s . This is mostly for his- 
torical comparison to [7j , since the current constraint on 
constant equation of state is w = —0.97 ± 0.08 [3l| (us- 
ing only geometric data independent of the sound speed) 
and so the effects of sound speed are suppressed due to 

lltCRfO, 

2) Early dark energy with constant speed of 
sound (cEDE). In order to allow for a period where w 
is further from —1 and so the sound speed has more in- 
fluence, we also consider a model with varying equation 
of state but constant sound speed. We choose the phe- 
nomenological early dark energy model of [32| but allow 
c s to be a free (constant) parameter. At early times w 
approaches in this model and so the value of c s can 
have observational consequences. The model parameters 
are the fraction of dark energy density at early times fl e 



3 



0.5 I 




o.oooi o.ooi o.oi o.i l o 0.2 0.4 0.6 0.8 1 

k [h Mpc-'] a 



FIG. 1: The deviation of the power spectrum of the matter 
density perturbations (Newtonian gauge) from the c s — 1 
case is plotted vs. wavenumber k. Three regions - above the 
Hubble scale (small k), below the sound horizon (large fc), and 
the transition in between - can clearly be seen. The models 
have w = —0.8 (deviations will be smaller for w closer to 
— 1) and constant sound speed as labeled. For the c s — 0.1 
case, we also show the result (dashed curve) in terms of the 
gauge invariant variable D g as defined in [3y] (in that work 
$ is equal to minus our <f>). This illustrates that the low k 
behavior is strongly gauge dependent. 



(this approaches a constant), the equation of state today 
wo, and c s . We call this generalization the cEDE model. 
Here 



^de (a) 
w{a) 



de 



a (i 



,-3i«o> 



fide "I - fim^^^ 

1 (i In fi^e (a) 

3[1 — fide (a)] din a 



n„ (i 



-3wq ■> 



(8) 



(9) 



where the current dark energy density VL^e = 1 — fi m - In 
this model, c s = const. We show an example of w{a) in 
Figure H 

3) Barotropic ("aether") dark energy models. 

The third model we treat is a particular case of the 
barotropic class of dark energy, where there is an explicit 
relation determining the pressure as a function of energy 
density. Several physical models for the origin of dark 
energy fall in this class, and have attractive properties as 
discussed below. 

Ref. [33| showed that all such viable models could be 
written as a sum of an asymptotic constant energy den- 
sity poo (with Wqc = —1) and a barotropic fluid, or 
aether, with positive equation of state wae > 0. The 



FIG. 2: The equation of state (lower three curves) and sound 
speed (upper three curves) as a function of scale factor are 
illustrated for two models. The aether model takes s = 3 
(solid curves) or s = 1 (dashed curves) and wo = —0.99; 
the early dark energy density Sl e is determined from these 
parameters. Note that the cEDE model (dotted curves, also 
taking too = —0.99, and here setting c 3 — 0) is a close match 
to the aether model. 



sound speed is completely determined by wae and has 
the property that c 2 s < wae- Moreover, to admit an 
early matter dominated era, wae{u <C 1) -> 0, and hence 
Cg(a <C 1) — > 0. We adopt the form wae = fia s so 



Pde{a) 
Pae{o) 
w(a) 



Pae(o) 



-3 3S(l-a s )/s 

PAE^a e py " 



(10) 

(11) 

Pae{o) , . 
WAE{a) — — (-12) 



c s (a) = WAE(a) - 



+ pae(o) f J oo + Pae{o) 

s wae{o) 



3 1 + wae (a) 



(13) 



where = pde.o ~ Pae.q- There are two free parame- 
ters in addition to the dark energy density today: /? and 
Pae,o ~ one less than in the cEDE case (we will fix s = 3 
usually). Note that the effective early dark energy den- 
sity Q e ss {pae.q/ Pm,o) e 3/3 / s and the present equation of 
state is wo = — 1 + {pae,q/ Pde,o) (1 + 0). As discussed 
by [33|, the barotropic model strongly ameliorates the 
coincidence problem, motivating why w w — 1 today. 

Our three models thus span constant w and constant 
c s , varying w and constant c s , and varying w and vary- 
ing c s (but with c s determined by w). We illustrate their 
equation of state and sound speed behaviors in Figure [21 
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We expect a cEDE early dark energy model with c s = 
to show the greatest effect of sound speed on the observ- 
ables. Since cEDE can look so much like the barotropic 
model, in w(a) and more approximately in c s , we do not 
treat the barotropic model separately in the following sec- 
tions, but rather consider it as a motivation for cEDE. 
The barotropic model possesses the advantage of having 
Cs = at early times (and wo ~ — 1 at late times) be- 
ing determined by physics rather than being adopted as 
phenomenology 



IV. IMPACT ON COSMOLOGICAL 
OBSERVATIONS 

We now consider angular power spectra of cosmological 
observables that are sensitive to the speed of sound of 
dark energy, with the aim of comparing the predictions 
to current observations (so we do not here include higher 
order correlations, leaving for future work such signatures 
and their effect on constraining non-Gaussianity). 
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A. Angular Power Spectra 

The matter density fluctuations, potential fluctuations, 
and the radiation field arc influenced by the dark energy 
sound speed as discussed in Sec. [TXJ From these we can 
form, and measure, the angular auto- and cross-power 
spectra. We consider the CMB temperature anisotropy 
power spectrum, the galaxy (or other large scale struc- 
ture tracer) overdensity field in a rcdshift bin i, and their 
crosscorrelation, giving the power spectra Cf^ , where 
{Xy} = {TT,Tgi, gigj}. See the Appendix for a review 
of how the power spectra relate to the potential power 
spectrum. 

Figure [3] shows a typical temperature power spec- 
trum. The signal from the sound speed dependence en- 
ters through the ISW effect, which is also plotted sep- 
arately in the figure. The extra power from the ISW 
effect arises from the decay of the potential as the dark 
energy impacts matter domination at late times; in the 
concordance model the cosmological constant dark en- 
ergy causes a decay in the potentials of about 25% be- 
tween the matter dominated era and the present. While 
the decay arises from the change in the expansion history 
due to the dark energy equation of state, it can be ame- 
liorated by increased dark energy clustering if the dark 
energy sound speed is small. Figure 2] illustrates the in- 
fluence of the sound speed. 

The ISW effect can be measured (341443! and one might 
hope to constrain the sound speed in this way. However, 
since the effect occurs only on the largest angular scales, 
cosmic variance swamps the signal. This is demonstrated 
in the left panel of Fig. [4] for a cosmic variance limited 
experiment scanning 3/4 of the whole sky. The right 
panel explicitly displays the low signal-to-noise for each 
multipolc, with the difference between c s — and c s = 



FIG. 3: CMB temperature power spectrum for w = —0.8 and 
c 3 — 1, explicitly showing the contribution of the late-time 
(z < 10) ISW effect. 



1 only amounting to S/N = 1 when summed over all 
multipoles. 

For the galaxy or matter density fluctuations, the dark 
energy sound speed can have a larger effect. Note that 
the dark energy perturbations themselves remain small 
relative to the matter inhomogeneities, despite a low 
sound speed having a dramatic effect on the dark en- 
ergy clustering. Figure [5] shows that on superhorizon 
scales the level of dark energy power is (1 + w) 2 rela- 
tive to the dark matter power (because at superhorizon 
scales the perturbations remain adiabatic and the ratio 
$de/o~dm = 1 + w). On subhorizon scales, the ratio 
depends strongly on the dark energy sound speed. For 
c s = 0, the ratio is scale independent in the subhorizon 
regime: during matter domination, one can show analyt- 
ically that then 



Pde 
Pdm 



1 



l-3w 



(matter dominated) (14) 



but this ratio becomes smaller by roughly a factor of two 
by today. For a canonical sound speed c s = 1, the dark 
energy power is strongly suppressed relative to the dark 
matter power, with the ratio scaling as fc -4 . 

The matter power spectrum itself, however, is affected 
by the dark energy sound speed through the poten- 
tial perturbations induced by the dark energy inhomo- 
geneities. Figure [6] shows in the left panel the absolute 
comparison of the dark matter and dark energy power 
(in contrast to the relative difference between the two in 
Fig. [5]). On this log scale one cannot see the influence of 
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FIG. 4: Left panel: CMB temperature power spectrum for c a = 0, and its difference from the c s = 1 case, are plotted for 
w — —0.8, along with the cosmic variance. Right panel: The signal relative to the noise (here just cosmic variance) is low, 
with the total summed over all multipoles S/N ~ 1.0. Compensating the difference between the models by varying the other 
cosmological parameters would make the S/N even smaller. 
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FIG. 5: The ratio of the dark energy to dark matter density 
power spectra (Newtonian gauge) is plotted for various values 
of constant w and c a . Although c s = gives dramatically 
more power on subhorizon scales than c 3 — 1, the direct ratio 
of the dark energy power to the matter power is negligible. 



the dark energy sound speed on the dark matter power, 
so the right panel plots the deviation with respect to the 
c s = 1 case. We see that the deviation due to c s = is 
at the percent level in the matter density power and the 
tens of percent level in the potential perturbation power. 

The density and potential are related through the Pois- 
son equation. For example, for w = —0.8 and c s = 0, the 
amplitude of the dark energy perturbations is about 4% 
of the dark matter perturbation (i.e. the power ratio is 
about 1.6x 10~ 3 on subhorizon scales as seen from Fig. [5]). 
According to the Poisson equation, Eq. (J7J, this trans- 
lates into about a 12% increase in <p g° m g from c s = 1 
to c s = 0, because today py>e ~ 3pm and because in the 
c s = 1 case the dark energy contribution to the Pois- 
son equation is negligible. Hence, as shown in the right 
panel of Fig [51 we get about a 25% increase in the power 
spectrum of 4>. 

Note that the (late) ISW effect is proportional to the 
change in potential A0 between matter domination and 
today. In the standard case, this decay is about 1/4 of the 
potential during matter domination and thus about 1/3 
of the potential today, i.e. Aip = <po — <j>MD ~ —\4>md ~ 
— h<j>o- Hence, the change in the potential at present of 
12% due to enhanced dark energy clustering corresponds 
to a change in the ISW effect of approximately 3x12% = 
36% (i.e. in [A0(c s = 0) - A0(c s = 1)]/A0(c s = 1)). 
This enhancement gives the ISW effect extra sensitivity 
to dark energy clustering relative to other probes. 

The matter density perturbation is of course also af- 
fected, but with only about a 1% increase in its ampli- 
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FIG. 6: Left panel: dark energy (lower four, thin curves) and dark matter (upper, thick curves) density power spectra for 
different choices of the dark energy equation of state and sound speed. Right panel: relative differences in the potential (<f>) and 
matter density (<5 m ) power spectra between c s — and c s — 1 (matter and dark energy perturbations in Newtonian gauge). 



tudc. This effect on the potential today through the Pois- 
son equation is therefore sub-dominant to the direct effect 
of the dark energy perturbation itself. 

Now that we have seen the basic effects of the dark 
energy sound speed and equation of state on the observ- 
ables, we consider the specific instances of the constant 
us model and cEDE model. We can already guess that to 
obtain reasonable constraint on the sound speed we will 
want a model that has as large a 1 + w and as small a c s 
as is consistent with the observations, for a substantial 
part of cosmic history. 

B. Estimating constraints in constant w Model 

We begin by estimating the chances of constraining the 
sound speed using the x 2 between two extremes: c s = 
and c s = 1. Since we consider angular power spectra 
and crosscorrelations of observables on the sky (labeled 
by capital letters below), \ 2 is in general given by 

X 2 = E E AC ? Y ( C ^)xVzw *C? W , (15) 
e {xy},{zw} 

where ACf Y is the difference in spectra between the two 
cases and the covariance is given by 

(Cov e ) XY ,zw = {2£ + 1)Uy ( G ? Z CT W + Gf w Gj z ) , 

(16) 

with 

C? Y =C? Y + N? Y , (17) 



where / s k y is the fraction of the sky that is observed, 
Cf- are the fiducial spectra and N* Y are the noise 
power spectra so that Ci are the observed power spec- 
tra that include the noise. (See the Appendix for further 
details.) For the \ 2 estimates of this section we only con- 
sider the CMB temperature power spectrum and we will 
consider the cosmic variance dominated limit where the 
noise power spectrum is much smaller than the fiducial 
power spectrum, Nf T = 0. Hence, Eq. (|T5|) simplifies to 

1 //\r TT \ 2 

x 2 = i/ sky E(^+ 1 )(^) ■ ( 18 ) 

Assuming Gaussian likelihood, the quantity \ 2 is 
equivalent to the signal to noise squared with which we 
can distinguish c s = 1 from our fiducial c a = if all the 
other parameters were known exactly. Since in reality 
we should marginalize over the other parameters as well, 
X 2 is an upper bound on the signal to noise squared for 
distinguishing the two sound speeds. Therefore if we find 
a low x 2 then there is little hope of constraining c s with 
the assumed dataset. To amplify the chances of detec- 
tion, we examine w = —0.8, since in the limit w — > — 1 
dark energy perturbations become irrelevant regardless 
of the value of the value of the sound speed; given that 
w = —0.8 is already an unlikely value given current data, 
the calculated signal to noise squared (S/N) 2 could be 
an overoptimistic estimate of the true value. 

Figure [5] confirms that the discrimination between 
sound speeds through the CMB temperature autocorre- 
lation is poor, as discussed in the previous subsection. 
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FIG. 7: Left panel: CMB temperature spectra for the early dark energy cEDE model with Sl e = 0.03, u>o = —0.8 is plotted for 
c s — and 1. The effect of changing the sound speed on the late ISW effect is a little stronger than in the case of ordinary 
w = —0.8 dark energy (also shown), but the major difference comes from higher £, where the early dark energy exhibits 
significant differences between c s = and c s = 1, while ordinary dark energy does not. Right panel: Signal to noise squared per 
mode for distinguishing c s = 1 from the c 3 — fiducial is plotted vs. multipole. The late ISW (treated as I < 21) contributes 
only (S/N) 2 = 1.8; including higher £, say all £ < 2000, gives (S/N) 2 = 8.8 x 10 3 . However, the differences at high £ can at 
least partly be compensated by varying other cosmological parameters. 



Cosmic variance swamps the difference between even the 
extremes, c s = and c s = 1, and the total (S/N) 2 w 1. 
Note that this took cosmic variance to be calculated from 
the most optimistic case, c s = 0, where the noise is sig- 
nificantly lower, so one truly cannot determine c s with 
the CMB temperature anisotropy despite all the most 
optimistic assumptions. 

The overall significance of the mere existence of the 
ISW effect (i.e. the \ 2 between the CMB power with 
the ISW effect artificially removed and the true CMB) is 
only {S/N) 2 ISW = 3.7. The potential decay in a model 
with dark energy sound speed c s = is a little less than 
half the contribution in the c s = 1 case, thus explaining 
the (S/N) 2 Aca=1 w (1/4) (S/N) 2 ISW = 1.0 quoted above. 
Thus the ISW signal in the CMB temperature spectrum 
is too blunt a tool to explore dark energy sound speed. 



C. Estimating constraints in cEDE model 



In the early dark energy case, we find that the ISW 
signal in both the CMB temperature autocorrelation and 
temperature-galaxy crosscorrelation is comparable to the 
signal in the case of ordinary dark energy (which typically 
has an energy density fraction relative to matter of ~ 
1CP 9 at CMB last scattering). However, there is another 
source of distinction. Dark energy in the cEDE model 
has w « at CMB last scattering; if in addition c s = 0, 
then cEDE behaves at early times just like dark matter, 
with significant clustering of the dark energy. This will 
affect not only the large scale, late time ISW contribution 
to the CMB but also the early Sachs- Wolfe effect and the 
acoustic peaks. 



We must go beyond the CMB temperature spectrum 
to consider the galaxy-galaxy power and temperature- 
galaxy crosscorrelation data. Rather than proceeding 
further with halfway measures such as calculating the 
signal to noise to determine whether we would be able to 
place to constraints on c s while fixing all other param- 
eters, we instead carry out a full likelihood analysis in 
Sec. H 



Therefore we expect a clearer observational signature 
of the sound speed than for ordinary dark energy. Fig- 
ure [7] shows the effect of changing the sound speed in 
the cEDE model. The CMB temperature autocorrela- 
tion alone delivers (S/N) 2 « 9 x 10 3 (for £ max = 2000). 
This seems more promising for constraining the sound 
speed, and again we proceed to a full likelihood analysis. 
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V. MEASURING THE SPEED OF DARKNESS 



To obtain accurate constraints on the dark energy 
sound speed we perform a Markov Chain Monte Carlo 
(MCMC) likelihood analysis over the set of parameters 
{\ogc s ,p dalk ,uj b ,uj c ,fl de ,T,A s ,n s }, where p dark is either 
w, in the constant w case, or {wq, f2 e }, in the cEDE case, 
u!b = £lbh 2 is the present physical baryonic energy den- 
sity density, tu c = tt c h 2 is the present physical cold dark 
matter energy density, fide is the present relative energy 
density in the dark energy, r is the reionization optical 
depth, A s the amplitude of primordial scalar perturba- 
tions (defined relative to a pivot scale of k = 0.05 Mpc -1 ) 
and n s is the spectral index of the primordial scalar per- 
turbations. Note that we choose log c s as the sound speed 
parameter because most of the sensitivity is at small val- 
ues of c s . 

For current data we include the CMB temperature 
power spectrum from WMAP5 44 1, the crosscorrclation 
of these temperature anisotropics with mass density trac- 
ers including the 2MASS (2-Micron All Sky Survey), 
SDSS LRG (Sloan Digital Sky Survey Luminous Red 
Galaxies), SDSS quasars, and NVSSJNRAO VLA All 

and the SDSS 
To break de- 



45 . 



Sky Survey) radio sources, following 
LRG autocorrelation function from 
generacies with background cosmology parameters and 
constrain the expansion history, we use the supernovae 
magnitude-redshift data from the Union2 compilation 

m. 

The MCMC package COSMOMC g| is used to cal- 
culate the joint and marginalized likelihoods. The re- 
sults for the marginalized ID probability distributions 
are shown in Fig. [5] for the constant equation of state 
case and in Fig. [5] for the early dark energy, cEDE case. 
Dotted lines show the distributions when one fixes c s = 1. 

In the constant w case, no constraint can be placed on 
the sound speed, as expected from our earlier arguments. 
In addition, the other parameter distributions are essen- 
tially unaffected by the value of c s . For the cEDE case, 
however, some preference appears for a low sound speed, 
c s < 0.1, and this propagates through to the other pa- 
rameters. Since early dark energy with a low sound speed 
acts like additional dark matter at early times, this allows 
a lower true matter density. 

It is intriguing to consider whether the apparent prefer- 
ence of current data for the ACDM model is merely a con- 
sequence of overly restricting the degrees of freedom of 
dark energy, and that instead a dark energy with dynam- 
ics (wq w —0.95), microphysics (c s rj 0.04), and long- 
time presence (f2 e « 0.02) could be the correct model. 

Figure [TO] shows the 68.3%, 95.4% and 99.7% confi- 
dence level contours in the w-log c s plane for the con- 
stant w model. We see that current data in this model 
prefer w ~ — 1 but are completely agnostic regarding c s . 
For the cEDE model, Fig. [TT] shows the joint probabil- 
ity contours among the dark energy parameters, in the 
u>o-logc s , f2 e -logc s , and fl e -wo planes, with all other pa- 
rameters marginalized. Here we sec that the model mcn- 




0.021 D.022 0.023 0.024 0.09 0.1 0.1 



1.03 1.035 1.04 1.045 




0.92 0.94 0.96 0.9S 1 



3 3.1 3.2 0.7 0.72 0.74 0.76 0.78 0.8 

iog[io'° a,] n. 




13.4 13.6 13.8 14 14.2 0.2 0.22 0.24 0.26 0.28 0.3 0.65 0.7 0.75 0.8 0.85 

Age/GYr fi m "b 




FIG. 8: Constant equation of state case, plotting the 
marginalized one dimensional probability distributions using 
data from supernovae (Union2), CMB (WMAP5), galaxy au- 
tocorrelation (SDSS LRG), and the cross correlation between 
large scale structure tracers (see text) and CMB temperature 
anisotropics. Solid lines are for the model with log(c s ) a free 
parameter (with a flat prior), whereas the dotted lines corre- 
spond to fixed c a — 1. 



tioned above, (u> , c s , Cl e ) = (—0.95,0.04,0.02), is com- 
pletely consistent with the data, as is the cosmological 
constant (—1,1,0). It will be interesting to see how the 
best fit evolves with future data. 



VI. CONCLUSIONS 

Current cosmological data are in excellent agreement 
with the standard ACDM universe, with equation of state 
w = — 1. Nevertheless, the current data are also consis- 
tent with a wide variety of richer physics. It is not clear 
that it is wise to assume that the physical explanation 
for dark energy in the universe is indeed given by re- 
striction to a spatially smooth, constant in time energy 
density: the cosmological constant. Even after allowing 
for dynamical dark energy, there could be further de- 
grees of freedom - "hidden variables" or microphysics - 
in the dark energy sector, harbingers of deeper physics 
that have not yet shown clear signatures in the data. An 
explicit search for these signatures, and thus the physics 
behind dark energy, should be near the top of the list of 
current efforts in cosmology. 

In this paper we search for degrees of freedom beyond 
quintessence by examining the influence of the sound 
speed of dark energy, and its resulting spatial clustering 
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FIG. 9: Early dark energy case, plotting the marginalized 
one dimensional probability distributions using data from 
supernovae (Union2), CMB (WMAP5), galaxy autocorrela- 
tion (SDSS LRG), and the cross correlation between large 
scale structure tracers (see text) and CMB temperature 
anisotropics. Solid lines are for the model with log(c s ) a free 
parameter (with a flat prior), whereas the dotted lines corre- 
spond to fixed c s = 1. 

of dark energy, on key observables and in current data. 
This extends earlier analyses, quantifying the effects on 
the dark matter and dark energy density perturbation 
power spectra, the potential power spectrum, and their 
crosscorrelation. Where possible, we give simple scalings 
with 1 + w and c s . We also explore models with time 
varying equation of state and sound speed. 

In the standard model with negligible dark energy at 
high redshift, the speed of sound is essentially not dis- 
tinguishable with current data (see Fig. [5]) because cur- 
rent data favor w ~ — 1, and the effects of clustering 
of dark energy vanish in this limit. As w gets further 
from —1, the influence of the sound speed increases; for 
models with w » at high redshift there is also the 
possibility of non-negligible amounts of early dark en- 
ergy density. Even just a couple percent of the total 
energy density in early dark energy dramatically im- 
proves the prospects for detecting dark energy cluster- 
ing. One can view the early dark energy fraction tt e as 
another degree of freedom to explore. Indeed, carrying 
out a MCMC analysis we find in Figs. [S] and [11] that 
a model with dynamics, microphysics, and persistence: 
(wo,c s ,Q e ) = (—0.95,0.04,0.02) is completely consistent 
with the current data (although A remains consistent as 
well). 

Discovery of the accelerating universe 12 years ago has 



FIG. 10: 68.3, 95.4 and 99.7% confidence level contours in 
the dark energy model with constant equation of state. The 
constraints are based on current data including CMB, su- 
pernovae, LRG power spectrum and crosscorrelation of CMB 
with matter tracers. 



propelled the physical interpretation of dark energy into 
one of the most important, exciting, and difficult prob- 
lems in physics. Although current observations indicate 
that the equation of state, as a constant or broadly aver- 
aged over time, is close to —1, this leaves considerable 
room for further physics, as demonstrated here using 
recent data. To go further we should explore all three 
frontiers of the dynamics w(a), the microphysics c s and 
spatial inhomogeneities, and the persistence fi e . 
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FIG. 11: 68.3%, 95.4% and 99.7% confidence level contours in the cEDE early dark energy model in the uio-logc s (left), fi e - 
logc s (middle) and Q, e -wo (right) planes. The constraints are based on current data including CMB, supernovae, LRG power 
spectrum, and crosscorrelation of CMB with matter tracers. 



Appendix A: Angular Power Spectra: Definitions 

Here we review how the observable power spectra 
of various quantities on the sky are related to the 
three-dimensional primordial power spectrum and the 
transfer functions. We consider the CMB temperature 
anisotropics and galaxy ovcrdcnsitics in redshift slices, 
or populations, labeled with the subscript j, and write 
the observables in direction n as line of sight projections 
along comoving radial coordinate x, 



X(n) = / rfx5 A (nx,r -x), 



(Al) 



with S (x, t) the "source term" as a function of comov- 
ing position and conformal time (to is the age of the 
universe in conformal time). Here X represents the ob- 
servable, which could be a galaxy overdensity gj in the 
jth redshift bin or a CMB temperature anisotropy T. For 
the galaxy overdensity gj , the source is 

S*@,t) = g(z(r)) " j( ^ T)) 6^ m (f,r), (A2) 
nf 

where rij(z)dz is the average angular galaxy density of 
galaxy population j in the redshift interval (z,z + dz), 
iVj- = J dzrij(z) is the total average angular galaxy den- 
sity of population j, and bj is the galaxy bias relative 
to the matter overdensity of bin j. The Hubble factor 
H(z) arises because the source was defined in terms of 
an integral over \ while rij(z)/n^ is normalized to unity 
in terms of an integral over z. 

For CMB temperature anisotropies, the (Fourier trans- 
form of the) source is given in Eq. (12) of [47|. The Inte- 
grated Sachs- Wolfe contribution to the CMB anisotropy 
is nonzero when the universe is not matter dominated, 
and thus the gravitational potentials <j) and ip are not 
constant. The ISW source is given by 



where dots denote derivatives with respect to conformal 
time. 

If we expand the anisotropy field in spherical harmon- 
ics, X(h) = J2e m a hn^rn{fi), the expansion coefficients 
are given by 

af m = J dnYUn)X(h) 

= (2n)-V 2 J dQY; m (n) J d 3 k J d X e^S x (k,r - x ) 

d 3 kY e * m (k) [ dxMkx)S x (k,T - X ), (A4) 



where we have Fourier expanded 



d 3 k 



s (f ' T) = j e S (fc > r) ' (A5) 

and we have used the Rayleigh plane- wave expansion 

e **ftx = A^i t h{k X )Yt m {k)Y lm {h) , (A6) 

where the je is the spherical Bessel function. We now 
write S x (k,r) — ^i(k) S x (k,r) where ipi(k) is the ini- 
tial potential perturbation and S x (k, r) is the source for 
ijji = 1, i.e. it is a transfer function. Due to the as- 
sumption of homogeneity, the transfer function does not 
depend on the direction of the wavenumber, but only on 
its magnitude k = \k\. The statistics of the initial per- 
turbations are given by 



^ bw (f,r)-0(f,r) + ^(x,T), 



(A3) 



(A7) 



where Pf (k) is the primordial potential power spectrum. 
Assuming statistical isotropy, the angular correlations 
between two quantities on the sky X and Y (where they 



11 



may or may not be the same) is given by the angular 
power spectrum 



( a fm a J'm') — Cf Y Su'S„ 

where, using Eq. (|A4I) . 



(A8) 



XY 



<PkP?(k) I dx3i(k X )S x (k,T - X ) x 



(2tt)2 

d x / Mkx')s Y (k,T -x') 



(A9) 



In this work, we are specifically interested in the com- 
binations {XY} = {TTjTg^gigj}, but Eq. (|A"9]) is the 
general expression for angular power or crosscorrelation 
spectra. 

When the sources S x and S Y vary slowly compared 
to the spherical Bessel functions in Eq. (|A9|) . the triple 
integral can to a good approximation be reduced to a 
single integral. Setting P(k) = P(k = (£ + 1/2) /x(z)) 
and using the asymptotic (for £ 3> 1) formula that 
(2/tt) / k 2 dku{kx)uik X ') = (1/x 2 ) S(x - X'), we find 



C, 



XY 



2n 2 



(£ + 1/2) 
x ft + 1/2 



d XX A} 



1/2 



S 



X 



.T 



) s y ( l -^,r im) 



where A(fc) = fc 3 P(fc)/(27r 2 ). We use the power spectra 
to calculate the X 2 (signal-to- noise) in Eq. (fT5|) . 



Finally, we need to specify formulae for noise in the ob- 
served spectra Cf Y . The covariances between the spec- 
tra are given by 



Cov(Cf r ,C, 



where 



tZW\ 



1 



(2£+ l)/ sky 



< r<xzr<YW 



XW^YZ^j 



(All) 



i~lX Y /~f 



XY 



N, 



XY 



(A12) 



Here / s k y is the sky coverage, Cf are the fiducial spec- 
tra and N XY are the noise power spectra. For the galaxy 
density fields, the white noise power spectra are given by 



nf 



and for the CMB it is given by 

Nj T = A T e £ ( f + 1 ) e FWHM/(81n2) 



(A13) 



(A14) 



where At is the sensitivity and #fwhm is the full width 
half max angle of the Gaussian beam. The noise cross 
power spectra can be assumed to vanish. 

The treatment of the covariances for actual data is 
typically more complicated than the above. In this paper, 
we use the covariances and treatment of the observables 
as given by the data packages in COSMOMC, 0, 0, 
|45| for the angular spectra, and the Union2 supernovae 
covariance matrix including systematics. 
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